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One percent of infants are born with congenital heart disease (CHD), which commonly 
involves outflow tract (OFT) defects. These infants often require complex surgeries, which 
are associated with long term adverse remodeling effects, and receive replacement 
valves with limited strength, biocompatibility, and growth capability. To address these 
problematic issues, researchers have carried out investigations in valve development and 
valve mechanics. A longstanding hypothesis is that flow-induced forces regulate fibrous 
valve development, however, the specific mechanisms behind this mechanotransduction 
remain unclear. The purpose of this study was to implement an in vitro system of outflow 
tract development to test the response of embryonic OFT tissues to fluid flow. A dynamic, 
three-dimensional bioreactor system was used to culture embryonic OFT tissue under 
different levels of flow as well as the absence of flow. In the absence of flow, OFT 
tissues took on a more primitive phenotype that is characteristic of early OFT cushion 
development where widely dispersed mesenchymal cells are surrounded by a sparse, 
disorganized extracellular matrix (ECM). Whereas OFT tissues subjected to physiologically 
matched flow formed compact mounds of cells, initated, fibrous ECM development, 
while prolonged supraphysiological flow resulted in abnormal tissue remodeling. This 
study indicates that both the timing and magnitude of flow alter cellular processes 
that determine if OFT precursor tissue undergoes normal or pathological development. 
Specifically, these experiments showed that flow-generated forces regulate the deposition 
and localization of fibrous ECM proteins, indicating that mechanosensitive signaling 
pathways are capable of driving pathological OFT development if flows are not ideal. 
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INTRODUCTION 

Outflow valve defects are involved in the most common 
Congenital Heart Disease (CHDs) (Schoen, 2008). Because valve 
defects are associated with stenosis and regurgitation may not 
become problematic until later in life, they are among the most 
costly and clinically relevant CHDs. For this reason, valve replace- 
ments, the second most common cardiac surgery (Combs and 
Yutzey, 2009), are being investigated to develop better therapies 
and treatments for the large population of both children and 
adults affected by valve disease. It is imperative that a portion 
of this research is focused on progressing a greater understand- 
ing of the developmental mechanisms that regulate healthy valve 
morphogenesis. 

It is well established that atrioventricular cushions (AVCs) and 
outflow cushions (OFCs) follow similar developmental mecha- 
nisms early in development (Hinton and Yutzey, 2011). In both 
cases, the valves are derived from the endocardial cushion tissue; 
however, OFCs lag slightly behind the AVCs and have a cellu- 
lar contribution from neural crest cells (NCCs) (Armstrong and 
Bischoff, 2004; Hinton and Yutzey, 201 1). In chick embryos, these 



NCCs undergo epithelial to mesenchymal transformation (EMT) 
around HH stage 1 1 before they begin migration into the dorsal 
OFT (Kirby and Hutson, 2010). By HH stage 21, these migrat- 
ing NCCs contribute to the OFCs (Webb et al, 2003) and follow 
similar signaling mechanisms as those found in AVCs, such as 
BMP signaling. Much like in the AVCs, there is evidence that 
flow-responsive signaling cascades play a role in the early develop- 
ment of OFCs. The misregulation of signaling pathways, such as 
KLF2/TGFf5 and RhoA/ROCK, lead to diseased semilunar valves 
that are characterized by increased, disorganized extracellular 
matrix (ECM) (Hinton et al., 2006). Further, studies have shown 
that reduced flow in the OFT results in impaired valve formation 
and septation (Hove et al., 2003), indicating that the connection 
between flow and development via specific signaling mechanisms 
is crucial to normal valve development. 

In a previous study, the impact of flow- induced forces on AVC 
development was investigated using an in vitro system of valve 
development with an extruded collagen tube scaffold (Tan et al., 
2012). Due to relatively large diameter of the original tube scaf- 
fold (4 mm), low velocity laminar flow profiles were required to 
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prevent the explant from washing away, restricting the range of 
potential flow magnitudes. For this reason, a new, molded (cast) 
collagen scaffold that contained a 0.65 mm restriction was used to 
achieve physiological levels of flow in an in vitro system of OFT 
development. In the study presented here, OFT explants were 
exposed to physiological and pathological levels of flow to inves- 
tigate the impact of flow-induced forces on the OFT cushion of 
the heart. Ultimately, this study contributes to the ever-increasing 
research efforts that seek to elucidate the developmental mecha- 
nisms that regulate healthy vs. pathological semilunar (SL) valve 
morphogenesis. 

MATERIALS AND METHODS 

CULTURE OF OFT EXPLANTS IN THE FLUID FLOW BI0REACT0R SYSTEM 

Avian embryos that have not yet reached hatching stage are 
not considered live animals. Nonetheless AVMA Guidelines for 
Euthanasia protocols were followed in carrying out these exper- 
iments. HH stage 25 OFT explants (approximately 10 per tube) 
were isolated by cutting the OFT at the heart/OFT boundary 
and then the OFT/aortic sac boundary. These tissues contain 
both conotruncal and intercalated cushions that participate in 
aorticopulmonary septation and OFT valve formation. The OFT 
explants were then butterflied open and seeded onto the nar- 
row channel section of the flow channel surface of the molded 
"hourglass" -like Coll scaffold (Figure 1). HH stage 25 OFT 
tissues were selected because NCCs have already migrated to 
the OFCs, and all cell types necessary for cushion develop- 
ment are present. The explants were allowed to attach to the 
scaffolds during a 4 day static culture period. The scaffolds 
containing OFT explants were then sutured onto the glass 
nozzles of the bioreactor. Reactors contained culture media: 
Dulbecco's Modified Eagle Medium (#50-013-PB, Corning 
Cellgro, Mediatech, Inc.) and 10%- Fetal Bovine Serum (FBS) 
(#F-0500-A, Atlas Biologicals, Inc.) with Penicillin Streptomycin 
(#P7794 and #S6501, Sigma- Aldrich), Amphotericin B (#A9528, 
Sigma-Aldrich), and Minocycline (#M9511, Sigma- Aldrich). A 
peristaltic pump (Ismatec Ecoline VC-MS/CA8-6, IDEX Health 
& Science) used with 1.85 mm inner diameter cassette tubing, 
was connected to the tubing circuit in the case of flow to gener- 
ate physiological or pathologically high levels of flow. "No-flow" 
control tubes were also sutured into bioreactors that were con- 
nected to a closed loop of stagnant media (not connected to the 
pump). 

AVCs and OFCs have similar dimensions and are exposed to 
similar blood flow rates early on as they follow common devel- 
opmental processes (Hinton et al, 2006; Liu et al., 2011). It is 
therefore likely that the level of shear stress is comparable between 
the two cushion types at early stages of development. A mean 
shear stress of 7.2 dyne/cm 2 was selected as the desired physio- 
logical condition based on previous studies of developing AVCs 
(Liu et al, 2009, 2011; Yalcin et al, 2011). To model the devel- 
opment of outflow tract stenosis, we chose a mean shear stress of 
16.6 dyne/cm 2 as a pathological flow condition. Since local flow 
profiles influence explant morphology we characterized the flow 
throughout the cast collagen tube by creating a three-dimensional 
computation fluid dynamics (CFD) model of the intraluminal 
culture media. The CFD model was then used to determine the 



flow rate needed to reach either 7.2 or 16.6 dyne/cm 2 in the 
seeding channel of the scaffold. 

Poiseuille's equation provided an analytical starting point for 
the CFD model; however, it was not sufficient to fully char- 
acterize the flow due to channel entrance length restrictions. 
The CFD model was implemented using COMSOL Multiphysics 
4.2a (Burlington, MA) to represent the molded collagen scaf- 
fold. The geometry was defined to match the dimensions of the 
mold's inner pins (scaffold lumen) and the structure was dis- 
cretized (Figure IB). A rigid wall assumption was applied and 
the media solution was considered Newtonian and incompress- 
ible with properties similar to that of water at 37 degrees C, 
(density, p, of 993.4 kg/m 3 and dynamic viscosity, fi, of 0.692 cP). 
The no-slip boundary condition was invoked at the wall and a 
pressure of 0 dyne/cm 2 was applied at the outlet. For the phys- 
iological condition, the inlet flow was adjusted until there was 
a mean shear stress of 7.2 dyne/cm 2 (flow rate of 24.3 mm 3 /s) 
or in the case of the pathological flow; a mean shear stress of 
16.6 dyne/cm 2 (flow rate of 48.7 mm 3 /s) was met at the center of 
the seeding channel. The Reynolds number at pathological con- 
ditions was found to be 137 indicating that the flow in this region 
was laminar for all flow conditions. 

The culture protocol used in these experiments is illustrated 
in Figure 1A. After the static 4-day seeding period, tubes were 
sutured onto bioreactors and the flow groups were exposed 
to introductory shear stresses of 3.4 dyne/cm 2 (flow rate of 
12.2 mm 3 /s) for 2 days as this initial rate was found be permissive 
for newly attached cushion tissues in previous experiments (Tan 
et al., 2012). Flow was then ramped up to match the mean physi- 
ological shear stress of 7.2 dyne/cm 2 (flow rate of 24.3 mm 3 /s) for 
2 days. To model the development of OFT stenosis, cultures were 
on reactors for an additional 2 during which the mean shear stress 
was increased to 16.6 dyne/cm 2 (flow rate of 48.7mm 3 /s). These 
pathological flow samples and 10 day no flow control cultures 
were then removed from reactors to represent the pathological 
case. 

3D RECONSTRUCTION OF SCAFFOLD WITH OFT EXPLANTS 

Cultured scaffolds with explants from all four groups (physio- 
logical flow with no flow control and pathological flow with no 
flow control, n = 2 for each group) were reconstruction three- 
dimensionally so that tissue morphogenesis in response to flow 
could be investigated. The samples were fixed in 4% PFA (in 
PBS buffer pH = 6.8) overnight, embedded in 5% agarose, and 
sectioned in 200 |xm increments using a vibratome (Oxford). 
Serial sections were mounted on slides, stained with HE, and 
imaged using a Zeiss Stemi DV4 Stereomicroscope (Carl Zeiss, 
Inc.). These serial images were then used to three-dimensionally 
reconstruct the tubes with explants in AMIRA (Visage Imaging, 
Inc.). After serial section alignment, the scaffold and explant tis- 
sue were labeled as separate materials (Figure 2, gold and purple, 
respectively). 

QUANTIFICATION TECHNIQUES 

To quantify the effects of flow vs. no flow on the OFT explants, 
qRT-PCR, confocal microscopy and RhoA activation ELISAs were 
performed. Real-time PCR primers were designed in Primer3 
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FIGURE 1 | Experimental design for flow-cultured outflow tract explants 
in molded scaffolds. (A) The total outflow tract (OFT) explant culture period 
(from time of dissection) was 8 days in the physiological case and 10 days in 
the pathological case. After the 4 day seeding period, a 2 day introductory 
flow period was used before flow was ramped to sub-physiological. After the 
2 day introductory period, some flow groups and no flow controls were 
removed and some remained on for an additional 2 days under physiological 



levels of flow. These final flow groups, and no flow controls, were then 
removed. (B) To determine the flow needed to achieve physiological or 
pathological levels wall shear stress at the center of the channel for outflow 
tract (OFT) explants (Purple), a 3-dimensional computational fluid dynamics 
(CFD) model was generated with precise geometry. Note the high fluid 
velocity in the seeding channel (red streamlines) and eddy formation 
downstream of the channel. 



(SimGene) and are listed in Table 1. The relative expression of 
genes of interest (GOI) was calculated [relative to an acidic ribo- 
somal binding protein (ARBP) standard] using the Pfaffl method, 
which normalizes the differences between the reaction kinetics 
of different primer sets (Pfaffl, 2001). In addition to qRT-PCR 
runs on the in vitro samples (n = 3 per group; 4 groups), tran- 
script levels were also quantified using freshly dissected OFTs 
from HH stages 25, 28, and 30 (n = 3 per group). These repre- 
sent the in vivo stages modeled over the time course of the in vitro 
experiments. 

Scaffolds were processed, sectioned, and prepared for confo- 
cal microscopy. Sections were stained for tenascin C (#AB19013, 
Millipore), periostin (#AB1404.1, Abeam), elastin (#AB21610, 
Abeam), Col6 (#5C6, Developmental Studies Hybridoma Bank), 
F-actin (rhoadamine phalloidin, R415, Invitrogen), fibronectin 
(#IST-9, Abeam), and P-Histone H3 (#9713S, Cell Signaling 
Technology, Inc., Boston, MA) using a concentration of 1:200 in 
BSA/PBS overnight at 4°C. "No primary antibody" controls were 



incubated in PBS only overnight at 4°C. The sections were then 
stained with Alexa-conjugated secondary antibodies (Invitrogen) 
and Dapi (#D21490, Invitrogen) overnight at concentrations of 
1:200 and 1:5000, respectively, overnight at 4°C. Sections were 
mounted and imaged using a Zeiss LSM 510 META confocal 
laser scanning microscope (Carl Zeiss). Images were collected 
using identical microscope settings per antibody so that relative 
expression could be determined between groups. Sample num- 
bers per group are listed in Table 2 for each antibody. Z-stacks 
were collected for each sample group using a Nyquist sampling 
approach and imported into the AMIRA software package. Once 
in AMIRA, the image stacks were rendered using identical thresh- 
old values for each channel across all samples. The module that 
yields voxels for each channel was then used to generate volumes 
for each protein of interest. This volume was taken as a ratio to the 
Dapi volume, which were to compare across all sample groups. In 
this way, the volume of protein is normalized by the number of 
nuclei in the tissue. 
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A Physiological Pathological 




B 

Physiological Pathological 




FIGURE 2 | Flow affects outflow tract explant morphology. (A) All 

outflow tract (OFT) sample groups were three-dimensionally reconstructed 
with labeled scaffold (gold) and explant tissue (purple). (B) Physiological 
flow groups exhibited a mound shape while pathological flow groups 
exhibited a smaller, leaflet morphology. In the case of no flow, tissue 
formed a loose network and occluded the lumen of the flow channel 
(comparisons are not to scale). Flow direction is indicated by the black 
arrow in three dimensional space. 



RhoA activity was quantified as described previously (Tan 
et al, 2012) using colorimetric ELISA assays (Total RhoA 
ELISA Kit, #BK150, Cytoskeleton Inc.; RhoA Activation Assay 
Biochem Kit, #BK124, Cytoskeleton Inc.). The percentage of 
active RhoA was subsequently determined using the dividend 
of active and total RhoA for each sample (physiological flow, 
pathological flow, and pathological no flow control, n = 3; 
physiological no flow control, n = 4). Activation assays were 
also performed on freshly dissected OFCs from HH stages 25, 
28, and 30 (n = 3 per stage) to quantify in vivo activation 
levels. 

TRANSMISSION ELECTRON MICROSCOPY 

Transmission Electron Microscopy (TEM) was performed to 
investigate the ultrastructure of the explants (n = 1 per group) 
were fixed in 2.5% glutaraldehyde in PBS (pH = 7.2) overnight 
at room temperature. Samples were subsequently washed in 
PBS, post-fixed in a 1% aqueous osmium tetroxide solu- 
tion, and dehydrated in a graded ethanol series. The sam- 
ples were then embedded in PolyBed 812 (Polysciences, Inc., 
Warrington, PA) and sectioned at 100 nm. Sample images 
were gathered with a JEOL 200CX TEM (JEOL USA, Inc., 
Peabody, MA). 



STATISTICAL ANALYSIS 

All data were expressed as Mean ± s.e.m. Statistical comparisons 
were performed using ANOVA or independent f-test by Prism 6 
(GraphPad Software, Inc., La Jolla, CA). P values <0.05 (*) were 
considered to be statistically significant. 

RESULTS 

FLOW-INDUCED FORCES AFFECT OFT EXPLANT MORPHOGENESIS 

The flow impacted the overall geometry of the OFT explant tis- 
sue (Figure 2). In the case of physiological flow, the explants 
formed a compact mound at the inflow side of the seeding 
channel. Explants exposed to pathologically high levels of flow 
and two additional days of culture formed a small, leaflet-like 
structure extending in the direction of flow from the seeding 
channel inlet and extending downstream through the channel of 
the scaffold. These pathological flow samples all had small areas 
of explant/scaffold interface. In no flow control groups, explants 
formed a loose, spider web-like network that bridged the entire 
luminal diameter of the scaffold channel. 

THE EFFECTS OF FLOW ON OFT EXPLANT FIBROUS ECM PROTEINS 

To investigate changes in transcript levels of fibrous ECM genes 
that are associated with cellular responses to mechanical stimuli, 
qRT-PCR was performed (Figure 3). Tenascin C, elastin, Coll, 
and Col6 transcripts were downregulated in flow cases com- 
pared to their no flow controls (Figures 3A,E,G,I, respectively), 
while periostin transcript was upregulated in response to flow 
(Figure 3C). Elastin, the only transcript that showed significant 
differences between physiological and pathological groups, was 
upregulated in the extended culture pathological groups (flow 
and no flow control) (Figure 3E). To enable the comparison of 
these in vitro results to changes transcript levels that occur during 
relevant stages of heart, qRT-PCR for freshly dissected OFCs from 
HH stage 25, 28, and 30 embryos. In vivo, all fibrous ECM protein 
transcripts appear to be upregulated between HH stage 28 and 
30 (Figures 3B,D,F,HJ). While transcriptional changes were seen 
in response to flow, to assess the changes in protein deposition, 
confocal microscopy was also performed. 

Confocal microscopic analysis indicated that both physio- 
logical and pathological flow promoted organized localization 
of several hallmark fibrous ECM proteins. Tenascin C and 
Col6 (Figures 4A,D, respectively) were localized on the inflow 
(upstream) side of the OFT cushions. Periostin (Figure 4B) was 
diffuse in the flow-cultured tissue with increased staining inten- 
sity along the luminal edge. Elastin (Figure 4C) staining was 
diffuse in the flow-cultured cushions with robust deposition 
on the outflow (downstream) side. All no flow control groups 
showed ECM staining of all four fibrous proteins throughout the 
tissue, with greater staining intensity at the scaffold-tissue inter- 
face, where some cells invaded the scaffold. This intrascaffold cell 
localization was not observed in any of the flow-cultured explants. 
Images were taken using the same microscope settings so that the 
mean fluorescence intensity (green and blue channels) could be 
determine and relative deposition quantified. 

Confocal image analysis indicated that flow, physiological 
or pathological, resulted in increased deposition of tenascin C, 
periostin, elastin, and Col6 when compared to no flow controls 
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Table 1 | Chick OFT experiment primer sequences. 



qRT-PCRGOl 


Forward 5'-to-3' 


Reverse 5'-to-3' 


ARBP 


ppt ptp atp tp a ptp tpp 
bb I b lb Alb I bA b I b Ibb 


APT — I — TP TPT PPP PTP TTA ATP 

Ab I I I b I b I bbb bib I I A Al b 


Tenascin C 


app ap a pat p pt ptp p a a pt 
Abb ALA LAb bb I b I b bAA b I 


TAP TPP PPP TP A PAP PTP AT 

I Ab I bb bbb I bA bAb blbAI 


Periostin 


PPA TPP TAT PAP ATP TP 

bbA I bb IAI bAb Abb Al b I b 


PPA AAP AAA PTP AATPAA PP 

bbA AAb AAA b I b AAI bAA bb 


Elastin 


ppp pap a ap pp a aat apt pt 
bbb bAb AAb bbA AAI Ab I b I 


PTP APP TTP PTP PPP AT — I — TA 

b I b Abb I I b b I b bbb Al I IA 


Col 1 a1 


Tpp tpt TP A TAP PAP PP A PT 

I bb I b I I bA lAb bAb bbA b I 


PTP TPP TPP PAP ATP APP TP 

b I b I bb I bb bAb Al b Abb I b 


Col6a1 


TPP APA TTP PPA APP ATT TP 

I bb AbA I I b bbA Abb All lb 


APP APP TTT TPP TTP PTP A A 

Abb Abb I I I I bb I I b b I b AA 


RhoA 


PAP PAP PPT PPA PTT PAP TA 

bAb bAb bb I bbA b I I bAb I A 


PPA TPP TPT PAP TPP APA A A 

bbA I bb lb I bAb I bb AbA AA 


Rac 


PPT — TTP AAA ATP TPP PTP PT 

bb I Mb AAA Alb I bb b I b b I 


PPT PPP TPT PTA APT PPT PA 

bb I bbb lb I b IA Ab I bb I bA 


FAK 


ATT PPT PPT APP A AP PTP PT 

Al I bb I bb I Abb AAb b I b b I 


APP PTP PPA APT — TP A TTP AP 

Abb b I b bbA Ab I I bA I I b Ab 


Vinculin 


ATP A A P PTP PAA AAPPAPPA 

Al b AAb L I b bAA AAb LAb bA 


T — 1 — T PPA PTP APP AAA TPP AP 

1 1 1 GLA b 1 b AbL AAA 1 LL Ab 


Paxillin 


PPP PTA A AT PPP APT PTP AT 

bLLLIAAAl bbLALI blbAI 


PT — T PTP PT — T PT — 1 TPP PAA PP 

CTT LTb L 1 1 b 1 1 1 bb bAA bL 


KLF2 


PPT TPT APP APA PA A APP PP 

bb I I L I ACL AbA LAA ALL Lb 


PAP PAP TPP PPP ATA APT PT 

LAb bAL 1 bb LLL Al A AL 1 b 1 


TGFpi 


GCT CTG TAC AAC AGC ACC CA 


ACG AAG AAG ATG CTG TGG CT 


TGFp2 


GAG AAA GCC AAC CAC AGA GC 


GGT ACA GCT CTA TCC GCT GC 


TGFp3 


TAA ATC CTC TCG TTG TGG GG 


CAC AAT GAG TTG GGC ATT TG 


CTGF 


CAG GTT TAG CGA CAG CAA CA 


GTA ACC TAA CTG CCG CAA GC 


aSMA 


GGG GAT GAT GCT CCA AGA GCA GTT T 


CCT CAG GTG CAA CAC GGA GCT CAT T 


a-Actinin 


CAC TTG CGG CAG TAT GAG AA 


CAG GGG TCA GAA TCT GGT TT 


Vimentin 


CCG ACA GGA TGT TGA CAA TG 


TCT TAG CAG CAA CGC TTT CA 



Table 2 | Sample number per group in OFT experiment confocal 
analysis. 



Antibody Physiological Phys. no Pathological Path, no 
flow flow control flow flow control 



Tenascin C 


n = 


4 


n - 


= 5 


n — 


5 


n = 


= 4 


Periostin 


n = 


7 


n - 


= 8 


n — 


8 


n = 


= 6 


Elastin 


n = 


8 


n - 


= 8 


n — 


4 


n - 


= 6 


Col6 


n = 


3 


n - 


= 3 


n — 


3 


n - 


= 3 


F-Actin 


n = 


3 


n - 


= 5 


n — 


7 


n - 


= 5 


Fibronectin 


n = 


8 


n = 


= 10 


n — 


6 


n = 


= 10 


P-Histone H3 


n — 


9 


n - 


= 3 


n — 


6 


n - 


= 5 



(Figure 5). Tenascin C levels were significantly higher in the 
physiological flow group than in the pathological flow group. 
Periostin deposition was significantly higher in the 8 day no flow 
group (pathological control) than in the 6 day no flow control. 
Elastin and Col6 exhibited similar expression patterns and nei- 
ther ECM protein showed a change between the physiological and 
pathological groups. 

THE EFFECTS OF FLOW ON CYTOSKELETAL ACTIN AND FIBRONECTIN 

While ECM proteins contribute to a cushion's mechanical prop- 
erties, their organization does not become pronounced until later 
in development when the stratified layers of ECM within the 
valves becomes pronounced (Hinton and Yutzey, 2011). Thus, 
at early stages in cushion development, the actin cytoskeleton 
and fibronectin (which acts as a ligand for integrin recep- 
tors i.e., cell/ECM binding), are likely to play a significant 
role in tissue stiffness. Accordingly, we analyzed F-actin and 
fibronectin deposition using confocal analysis. Physiological 



levels of flow increased F-actin and fibronectin deposition 
compared to the physiological no flow controls while pathological 
flow decreased F-actin and fibronectin protein expression com- 
pared to the pathological no flow controls (Figures 6A,B). F-actin 
and fibronectin were diffuse throughout the tissue with higher 
staining intensity at the scaffold interface in all sample groups 
(Figures 6C,D). 

CHANGES OFT EXPLANT RhoA SIGNALING 

Previously, we investigated PJioA as a possible mechano- 
transduction-signaling molecule in response to flow-induced 
forces. In these studies, we found that RhoA has a role in flow- 
induced ECM development in AVCs (Tan et al., 2012). Here, we 
found little to no difference in RhoA activation was seen between 
flow groups and their no flow control counterparts in OFT 
explants (Figure 7A). RhoA activation was higher in the phys- 
iological groups than the pathological groups. However, RhoA 
activity in freshly dissected OFTs was found to steadily decrease 
from HH stage 25 to 30 (Figure 7B). 

The small GTPases regulate actin reorganization and inter- 
act with a number of other pathways. So, the transcript lev- 
els of related genes were investigated (Figure 8). RhoA, FAK, 
and vinculin transcript levels were higher under physiological 
flow than the no flow controls (Figures 8A,E,G, respectively). 
Rac and vinculin transcript decreased with culture time in 
flow and no flow control groups (Figures 8C,G). The RhoA 
flow groups showed decreased transcript over time (Figure 8A) 
and the paxillin no flow control transcripts decreased with 
time (Figure 81). No changes in FAK or Rac transcripts were 
seen from HH stage 25 to 28 (Figures 8D,F, respectively) 
while RhoA, vinculin, and paxillin all decreased from HH 
25 to 28 and increased from HH 28 to 30 (Figures 8B,H,J, 
respectively). 
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FIGURE 3 | Flow affects fibrous extracellular matrix protein transcript 
levels in outflow tract. Quantitative Real Time PCR was performed on the 
four culture groups and on freshly dissected outflow cushions (OFCs) at 
three stages in development (HH stage 25, 28, and 30). (A,B) Tenascin C 
transcript level was higher in the absence of flow and increases between HH 
stage 28 and 30. (C,D) Periostin transcript level was higher in the presence of 
flow and appears to steadily increase from HH stage 25 to 30. (E,F) Elastin 



transcript level was higher in the absence of flow and in pathological groups 
and transcript level appears to increase between HH stage 28 and 30. (G,H) 
Type 1 collagen (Coll) transcript level was higher in the absence of flow with 
no difference between physiological and pathological groups. Coll transcript 
appears to increase at HH stage 30. (I,J) Type 6 collagen (Col6) followed a 
similar trend to Col 1 but with no statistical significance and Col6 transcript 
levels appear to steady increase from HH stage 25 to 30. CP < 0.05). 
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FIGURE 4 | Flow affects fibrous extracellular matrix protein localization 
in outflow tract. Outflow tract (OFT) explant samples were stained with 
Dapi (blue) or for several hallmark extracellular matrix (ECM) proteins 
(green): (A) Tenascin C (B) Periostin (C) Elastin (D) Type 6 Collagen (Col6). 
In all no flow controls, ECM staining was throughout the tissue network, 
with a greater staining concentration at the scaffold interface. In each case 
of flow, tenascin C, periostin, and Col6 were localized toward the inlet side 
of the cushion (A,B,D) while elastin staining presented throughout the 
cushions with greater concentration on the outlet side. Scale bars are all 
100 |im. The dotted lines indicate the scaffold wall. 



RhoA has also been connected to shear-responsive signaling 
molecules including KLF2 and the TGFf5 proteins. A downstream 
target of RhoA in this pathway is CTGF, which is associated with 
fibrous ECM protein production. TGF|3 1, 2, and 3 transcript lev- 
els were all higher in the presences of flow (physiological and 
pathological) while KLF2 transcripts were increased in only phys- 
iological flow samples (Figures 9A,C,E,G, respectively). In the 
physiological groups, CTGF transcription was downregulated in 
response to flow (Figure 91). In vivo OFTs were analyzed and 
showed that KLF2 and TGFfU transcripts decreased from HH 
stage 25 to 30 (Figures 6, 9 B,D, respectively) and CTGF increased 
from HH stage 25 to 30 (Figure 9J). No significant change in 
TGFf52 and 3 transcripts were seen between HH stage 25 and 30 
(Figures 9H,F, respectively). 

Finally, to look further at possible cushion cell phenotype, 
aSMA, a-actinin, and vimentin transcripts were analyzed. aSMA 
and a-actinin flow group transcript levels both decreased with 
time (Figures 10A,C, respectively). Vimentin transcript levels 
were upregulated in response to flow with no significant change 
between the two flow groups (i.e. time cultured under flow) 
(Figure 10E). aSMA transcript decreased from HH stage 25 to 
30, while a-actinin transcript increased from HH stage 25 to 30 
(Figures 10B,D, respectively). An increase in vimentin transcript 
was seen from HH stage 28 to 30 (Figure 10F). 

The normalized relative quantity (NRQ) ratios between dif- 
ferent group combinations are listed in Table 3 (in vitro culture 



groups) and Table 4 ("in vivo" samples). In each table, a color 
gradient is applied to compare the relative magnitudes of each 
NRQ ratio. Bright red is indicative of a significant decrease (P < 
0.05) while bright green is indicative of a significant increase 
(P < 0.05). White cells indicate little to no difference between the 
compared groups. 

FLOW ALTERS OFT CELL DEATH AND PROLIFERATION 

To look at the ultrastructure of the OFC culture groups, trans- 
mission electron microscopy (TEM) was carried out. Numerous 
sarcomeres with associated z-bands were observed (arrowheads 
Figure 11). These cells, which were present in all groups, also 
had centrally located nuclei indicating that they were cardiac 
myocytes (Figures 11 A-D). In contrast, endothelial cells were 
only observed in the physiological flow and no flow control 
groups (Figure 11A,B, respectively). In physiological flow group, 
membranous protrusions were observed in on the luminal side of 
these endothelial cells (Figure 11 A). However, no microtubules 
were associated with these protrusions in any of the groups 
(Figure 11 A). A loose network of ECM was observed in both 
flow groups (Figures 11A,C). All groups showed evidence of cell 
debris; however, this was more wide spread in the supraphysio- 
logical, pathological flow group (Figure 11C). 

Because p-Histone H3, is expressed in the myocardium during 
pathological tissue apoptosis (de Vlaming et al., 2012), sections 
from the four culture groups were stained with p-Histone H3 and 
Dapi and the number of positive nuclei vs. total nuclei per section 
was calculated. Pathologically high levels of flow (supraphys- 
iological) were associated with increased cellular proliferation 
compared to the no flow control and the physiological flow group 
(Figure HE). 

DISCUSSION 

The goal of the experiments presented in this study was to 
characterize OFC cellular response to predicted physiological 
or pathologically high (supraphysiological) flow-induced forces. 
This response was characterized as normal or abnormal based on 
tissue morphology, fibrous ECM protein expression, and local- 
ization, comparison of many different gene transcripts associated 
with the post-EMT developmental stage, and subcellular ultra- 
structure. Particular changes in all five of these areas are associated 
with valve diseases, the most common CHDs. As valve disease is 
problematic and costly to both children and adults almost equally 
(Schoen, 2008), understanding the mechanistic progression of 
valve pathogenesis is of extreme importance. 

These experiments found that HH stage 25 OFT explants cul- 
tured under predicted physiological flow lead to the formation 
of a compact mound, similar to the cushions that are seen during 
early stages of development. Pathologically high flow cultures lead 
to thin leaflet formation with a reduced scaffold interface area. 
This "leaflet detachment" may represent failure of the tissue's cel- 
lular architecture under pathologically high loading conditions. 
In the absence of flow, a loose network of cells formed that 
occluded the lumen. Cells in the no flow controls invaded the 
walls of the collagen scaffold, indicative of a possible EMT re- 
initiation. This was also observed in our previous report using 
AV explants (Tan et al, 2012). 
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FIGURE 5 | Flow affects fibrous extracellular matrix protein expression in 
outflow tract. Confocal microscopy was performed for outflow tract (OFT) 
explants at the same settings and the mean fluorescence intensity of each 
fibrous extracellular matrix (ECM) protein was analyzed relative to Dapi 
staining. Expression quantifications were averaged for each of the four 
culture groups. (A) Tenascin C was upregulated under physiological flow 



compared to the no flow control and pathological flow. (B) Periostin was 
upregulated under physiological flow compared to the no flow control. The 
extended no flow culture period also appeared to upregulate periostin. (C) 
Elastin was upregulated under physiological flow compared to the no flow 
control and no differences were seen between physiological and pathological 
groups. (D) Col6 followed a similar trend to elastin. (*P < 0.05). 



A interesting difference between these studies with OFT tis- 
sues and our previous AV cushion studies are that fibrous ECM 
gene transcripts were upregulated by flow in the AV cushions, 
but not in OFT cushions. However, fibrous ECM protein depo- 
sition, which increases during normal development, was found 
to increase in response to both levels of flow in OFT tissues. 
The no flow controls, by comparison, had sparse deposition of 
these proteins, which is, again similar to our AV cushion stud- 
ies. This indicates that flow plays a central role in stimulating the 
fibrous ECM protein deposition the OFT post-transcriptionally. 
Collagen and tenascin C have been found to colocalize early 
on in OFC development (Tan et al, 2011), and in these exper- 
iments, Coll and tenascin C colocalized only in response to 
flow. Elastin also increased in response to flow and accumulated 
throughout the OFT explant cultures, with higher concentration 
on the downstream side. Simulation results of flow through the 
molded collagen scaffold indicated downstream eddy formation 
(Figure IB), so this elastin localization may be associated with 
localized stresses due to the eddies. 

Previously, we found increased fibrous ECM protein expres- 
sion due to extremely low levels of laminar flow past AV explants 
was attributed to increased tissue stiffness (Tan et al., 2012). While 
ECM stratification is extremely important to valve mechanics late 
in development, fibrous ECM proteins such as collagen are rel- 
atively sparse early in development and the cellular cytoskeleton 
becomes important in determining overall tissue mechanics. F- 
actin and fibronectin were found to localize throughout the 



explant tissue in the 4 culture groups, and exhibited higher 
staining concentration at the scaffold wall interface. Both F-actin 
and fibronectin expression were higher under physiological flow, 
and pathological flow appeared to reduce expression. This data 
indicates that physiological flow increases actin polymerization 
and thus may promote the mechanical contribution of the 
cytoskeleton, while pathologically high flow has the opposite 
effect. 

Though flow increased deposition of fibrous ECM proteins in 
OFT explants, flow decreased tenascin C, elastin, Coll, and Col6 
transcripts and increased the periostin transcript. The expression 
of fibrous ECM protein such as collagen is complex with numer- 
ous post-transcriptional regulatory mechanisms. Importantly in 
our experiments, flow was found to increase protein deposition 
in the OFT, which is critical in in vivo development. In our previ- 
ous report using AV explants, these transcripts were all increased 
(Tan et al., 2012). This indicates that flow regulation of these tran- 
script levels may be more complex in the OFT. Perhaps different 
waveforms are necessary to drive increases in transcript levels but 
the steady state flows delivered in our current system are capable 
of driving protein deposition. 

Flow stimulated regulation of periostin appears to be different 
than the other fibrous ECM proteins investigated here. Periostin, 
which is necessary for collagen fibrilogenesis, is responsive to 
TGFf53, a gene that exhibits increased transcript in response 
to flow (Norris et al, 2007, 2009; Hinton and Yutzey, 2011). 
Previous in vitro studies using AV explant cultures showed that 
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FIGURE 6 | Flow affects cytoskeletal actin and extracellular fibronectin in 
OFT cells. Mean fluorescence intensity relative to Dapi was quantified from 
confocal imaging under the same setting for each outflow tract (OFT) explant 
culture group. (A) F-actin expression was upregulated under physiological 
flow compared to no flow controls. (B) Fibronectin expression was 



upregulated under physiological flow compared to no flow controls and 
extended no flow culture time appeared to upregulate expression. (C) F-Actin 
and (D) fibronectin expression both exhibited diffuse staining in all sample 
groups with increased concentration at the scaffold-tissue interface. Scale 
bars are all 100 \im. The dotted lines indicate the scaffold wall. (*P < 0.05). 




FIGURE 7 | Flow alters RhoA activation in outflow tract explants. Total 
RhoA and active (GTP-bound) RhoA ELISAs were performed to quantify the 
ratio of active vs. total RhoA. (A) In the physiological group, more RhoA 



appeared to be activated in the absence of flow (no flow control), and RhoA 
activation was decreased in the pathological groups. (B) RhoA activity 
appears to decrease steadily from HH stage 25 to 30. CP < 0.05). 



periostin was not induced by the RhoA pathway, and it was 
proposed that regulation of periostin was regulated by other path- 
ways (Tan et al., 2012). Periostin was; however, upregulated by 
LPA, so possibly, another Rho family GTPase, such as Rac or CDC 
42, is involved in periostin regulation. All fibrous ECM protein 
transcripts were found to increase in OFTs from HH stage 28 to 
30. Transcripts of the flow-responsive genes, KLF2, TGFfil23, and 
RhoA increased in response to physiological flow, while CTGF 



decreased in response to flow. This opposing trend is the same 
in the in vivo cases (i.e. if CTGF transcripts are upregulated, 
then KLF2, TGF|3123, and RhoA transcripts are downregulated). 
While RhoA/actin dynamics are necessary for CTGF and tenascin 
C induction (Asparuhova et al., 2009), they are also necessary 
for other cell processes such as cytoskeletal actin organization. 
Possibly, at this early stage in development, fibrous ECM pro- 
teins, while they are characteristically more organized, are not 
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FIGURE 8 | Flow affects adhesion signaling transcript levels in outflow 
tract. Quantitative Real Time PCR was performed on the four culture groups 
and on freshly dissected outflow cushions (OFCs) at three stages in 
development (HH stage 25, 28, and 30). (A,B) RhoA transcript level was 
higher under physiological levels of flow and appears to decrease from HH 
stage 25 to 28 and increase from HH stage 28 to 30. (CD) Rac transcript 
level was higher in the physiological sample groups (flow and no flow control) 
and no significant differences were seen from HH stage 25 to 30. (E,F) FAK 



transcript level was higher in the case of flow and no significant differences 
were seen from HH stage 25 to 30. (G,H) Vinculin transcript level was higher 
for each flow case compared to no flow controls and transcript was greater in 
the physiological groups compared to the pathological groups. Transcript 
appears to, like RhoA, decrease from HH stage 25 to 28 and increase from 
HH stage 28 to 30. (I,J) Paxillin transcript level decreased with no flow 
culture time and appears to follow the same trend as RhoA and vinculin from 
HH stage 25 to 30. (*P<0.05). 
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FIGURE 9 | Flow affects shear-responsive signaling transcript levels in 
outflow tract. Quantitative Real Time PCR was performed on the four 
culture groups and on freshly dissected outflow cushions (OFCs) at three 
stages in development (HH stage 25, 28, and 30). (A,B) KLF2 transcript 
level was higher in the absence of flow and increased with increased 
static culture time. Transcript appears to decrease between HH stage 25 
and 28. (C,D) TGF(S1 transcript level was higher in the presence of 
physiological flow compared to the no flow control. TGFf51 transcript 



steadily decreases from HH stage 25 to 30. (E,F) TGF(52 transcript level 
was higher in the presence of flow and appears to decrease from HH 
stage 25 to 28 and then increase from HH stage 28 and 30. (G,H) 
TGF(S3 transcript level was higher in the presence of flow with no 
difference between physiological and pathological groups. Transcript 
amount does not appear to change from HH stage 25 to 30. (1,J) CTGF 
transcript level was highest in the physiological no flow control group 
and appears to increase from HH stage 28 to 30. i*P < 0.05). 
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FIGURE 10 | Flow affects cytoskeleton remodeling transcript levels in 
outflow tract. Quantitative Real Time PCR was performed on the four culture 
groups and on freshly dissected outflow cushions (OFCs) at three stages in 
development (HH stage 25, 28, and 30). (A,B) aSMA transcript level was 
highest in the presence of physiological flow and appears to decrease between 



HH stage 28 and 30. (C,D) a-Actinin transcript level was the same between 
flow groups and their controls, but transcript was greater in the physiological 
groups. Transcript appears decrease from HH stage 25 to 28 and increase from 
HH stage 28 to 30. (E,F) Vimentin transcript level was higher in the presence of 
flow and appears to increase from H H stage 28 to 30. (* P < 0.05). 



the primary cellular output in response to external mechanical 
stimuli. 

Transcription of migration-associated molecules including 
FAK, RhoA, vinculin, and otSMA (a myofibroblast marker) 
increased when explants were exposed to physiological flow. 
However, when cultured under pathological conditions, RhoA 
activity significantly decreased and pathological flow decreased 
aSMA transcription. FAK is thought to regulate polarized RhoA 
and Rac activity during cell migration, so increases in migratory 
molecule transcription may indicate a healthy remodeling pro- 
cess in response to flow-induced loading (Li et al., 2005). Further, 
RhoA, Rac, and FAK transcripts do not change from HH stage 
25 to 30, indicating that these cytoskeletal remodeling processes 
are ongoing and necessary throughout development. TEMs also 
showed some possible signs of cell migration in the physiolog- 
ical flow group, including membranous protrusions associated 
with cushion development (Markwald et al., 1975). Pathological 



levels of flow appeared to drive cell apoptosis, indicated by qual- 
itative assessment of dead cell debris in TEMS, and quantitative 
increased proliferation, identified by an increase in p-Histone H3 
positive cells, which are both characteristic of valve pathology 
(Pfaffl, 2001). 

The finding of sarcomere containing cells in the TEM stud- 
ies indicates that myocardialization of the OFT explants was 
in progress. Cardiac myocytes were found in all study groups, 
but appeared to be more abundant in the flow groups. Thus, 
flow may have a role in this step of OFT development. The 
endothelial cell protrusions found on the luminal surface of 
the physiological flow cultures provide a potential mechanism 
of flow sensing in these cells. We looked for primary cilia in 
these structures as they have been reported by others (Egorova 
et al., 2011). However, none were found in our studies. This 
may be due to the levels of flow used in this study, as these 
cilia have been reported to retract in regions of high flow. 
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Table 3 | Comparison of NRQ ratios in in vitro oft experiments. 

qRT-PCR Physiological Pathological Physiological Physiological vs. 
GOI flow vs. flow vs. vs. pathological pathological 

no flow no flow flow no flow 



Tenascin C 




Periostin 


2.07 


Elastin 


0.38 


Coll 


0.20 


C0I6 


0.50 


RhoA 


2.29 


Rac 


1.04 


FAK 


2.02 


Vinculin 


1.24 


Paxillin 


0.76 


KLF2 


1.31 


TGFP1 


1.90 


TGF02 


1.22 


TGFP3 


1.49 


CTGF 


0.61 


ctSMA 


1.12 


a-Actinin 


1.00 


Vimentin 


1.42 



0.44 

2.72 
0.41 
0.21 

0.40 

0.95 
1.13 
1.53 
1.30 
0.97 
0.90 
0.93 
1.95 
1.30 
1.19 

1.05 
1.97 



0.69 
0.90 
0.39 
0.86 
1.05 
2.48 
1.33 
0.87 
1.33 
1.19 
0.92 
0.83 
0.82 
0.84 
0.74 
6.36 
1.91 
1.04 



1.41 
1.1E 
0/ 
0.91 
0.83 
1.03 
1.45 
0.66 
1.39 
1.52 
0.64 
0.41 
1.31 
0.73 
1.45 
1.45 
2.00 
1.44 



Table 4 | Comparison of NRQ ratios in OFT for three OFT 
developmental stages. 



qRT-PCR 


Stage 28 


Stage 30 


Stage 30 
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vs. stage 25 


vs. stage 26 


vs. stage 25 
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Clearly a more in depth study is required to clarify these 
questions. 

The study presented here illuminates the importance of flow- 
induced loading on cushion developmental mechanisms that 
ultimately regulate healthy vs. pathological valve mechanics. 
Numerous CHDs are attributed to OFT malformation, involv- 
ing both the valves and the arterial poles (Lucitti et al, 2005, 
2006; Schoen, 2008; Freund et al, 2012). Using an in vitro system 
of OFT explant development, physiological levels of flow were 
shown to promote cellular processes characteristic of normal 




' ' ' 



4 



Physiological Flow 200 micron. n 0 plow (Phys. Control) 





200 microns 



Pathological Flow 




No Flow (Path. Control) 



Proliferation 



0.4 



« 0.3 



2 0.2 



0.1 




Physiological 



Pathological 



FIGURE 11 I Flow alters outflow tract explant cell death and 
proliferation. Transmission Electron Microscopy was performed to 
investigate subcellular architecture. All samples contained myocyte cells as 
indicated by the presence of sarcomeres (arrowheads) and centrally located 
nuclei (A) Samples subjected to physiological flow exhibited an endothelial 
layer lining the flow side of the cushion, as well as ECM deposition and 
membranous protrusions (arrow). (B) Physiological no flow controls had 
some dead cell debris beneath the luminal cell layer. (C) Pathological flow 
groups did not show a luminal endothelial layer but showed MCs 
surrounded by ECM and cell debris. (D) Pathological no flow controls were 
similar to physiological no flow controls. Scale bars are all 200 |xm. (E) 
Confocal microscopy was performed on sample groups for P-Histone H3 as 
a proliferative marker. P-Histone H3 (+) nuclei relative to total nuclei (Dapi) 
were quantified and indicated increase proliferation in the case of 
pathological flow. (*P < 0.05). 



valve development, while supraphysiological levels of flow were 
shown to impede tissue development and drive pathobiology. 
The data presented here provides a necessary contribution to the 
growing field of heart valve engineering. By shedding light on 
the role of flow-induced forces on these important developmental 
processes, new possibilities can be laid out in the development of 
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better therapies and replacements for both CHD and non-CHD 
attributed valve disease. 

LIMITATIONS OF THE STUDY 

The central goal of these studies was to determine the role that 
fluid flow plays in the fibrous development of the aorticopul- 
monary septal and valvular tissues. A significant limitation of our 
study is that we were only able to deliver steady flow in this sys- 
tem. In vivo flow dynamics are much more complicated during 
morphogenesis of the outflow tract. Vortical, pulsatile, and even 
regurgitant flows are present during development of the arte- 
rial pole of the heart. To better model these fluid dynamics, we 
are developing flow systems that can deliver these more com- 
plex flows to cardiovascular tissues as they are likely to influence 
morphogenesis. 

Another limitation of this study is in the experimental design. 
We chose to model the development of stenosis after an initial 
"normal phase" of development. Thus, pathological studies lasted 
two days longer than the physiological studies. This means that 
the differences that we found between the supraphysiological and 
physiological conditions could be due to the lengthened culture 
time and not the increase in shear stress. We plan to use our more 
complex flow bioreactor system to parse out these questions and 
many others with greater resolution and precision. 
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